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Abstract 

The quasi-stationary nonequilibrium distribution function of an independent electron gas inter- 
acting with a medium, which is at local thermal equilibrium, can be obtained by entropy production 
rate minimization, subject to constraints of fixed moments. The result is not restricted to the region 
near equilibrium (linear response) and provides a closure of the associated generalized hydrody- 
namic equations of the electron gas for an arbitrary number of moments. Besides an access to 
far from equilibrium states, the approach provides a useful description of semi-classical transport 
in mesoscopic conductors, particularly because macroscopic contacts can be naturally taken into 
account. 

PACS numbers: 05.60.-k, 71.10.Ca, 73.23.Ad 
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Introduction - Electron transport in matter is often described by the Boltzmann transport 
equation (BTE) for the distribution function /, from which the (semi-classical) transport 
properties can be calculated Kohler 2] proved that the stationary solution of the lin- 
earized BTE satisfies a variational principle for the entropyproduction rate, which has been 
widely used to determine linear transport coefficients jl, ^, 4, 5]. Schlup and Jones 
gave arguments for the validity of Kohler's principle beyond linear response. Note that a 
linear BTE does in general not imply a restriction to the near-equilibrium, i.e., linear re- 
sponse, region. First, provided the linearity is not due to linearization of a nonlinear BTE, 
the linear BTE may be valid for large deviations from equilibrium. Secondly, because forces 
appear as coefficients of a term linear in /, the resulting currents are generally nonlinear 
functions of the forces. Below, a method is introduced that provides / as a function of its 
moments by minimization of the total entropy production rate. The result is not restricted 
to the near-equilibrium case (i.e., to linear response), and serves as a closure of generalized 
hydrodynamic equations for the moments. Recently, it has been shown that an analogous 
method applied to a photon gas in matter at local thermal equilibrium, where the BTE is 
exactly linear, describes nonequilibrium radiation satisfactorily well also far from equilib- 
rium [8|. 

After introducing the method, a few examples will be discussed, including metallic and 
nonmetallic electric conduction, and low-frequency transport in mesoscopic conductors. 



Basic formulation of the problem - Consider electrons in d{= 1, 2, 3)-dimensional space, 
with space and velocity vectors x and v, respectively. They may interact with a medium 
at local thermal equilibrium, such that the electron distribution function /(x, v) obeys the 
linear BTE 

+ V,/ + a- V,/ = /:(/o-/) . (1) 

The force m a = e V^^f/ acting on the particles with effective mass m and electron charge — e, 
is related to the gradient of the electric potential f/(x), which may be determined later (on 
the hydrodynamic level) self-consistently from the Poisson equation. The positive definite 
and self-adjoint linear (integral-) operator C describes the interaction of the electrons with 
the medium. The terms —Cf and Cfo can be interpreted, respectively, as absorption of 
electrons and emission of electrons equilibrated with respect to the medium (e.g., traps, 
contacts in mesoscopic systems, etc.), which has local temperature T and (electro-) chemical 
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potential jj,. The index indicates the Fermi equihbrium distribution /o = {exp(-^) + l}~^, 
where H{v) is the particle energy, f =| v |, and k is the^Boltzmann constant. Equation ([1]) 
applies to a large class of systems in solid state physics [1], presumes micro-reversibility, and 
is usually solved within linear response, i.e., in first order of 5/ = / — /q. Typical standard 
methods are the BGK approximation [9], or Kohler's variational approach with the help of 
trial functions. However, if the linearity is not due to a linearization of a nonlinear scattering 
term but remains valid for larger deviations of / from /o, a nonlinear (far from equilibrium) 
solution is apposite. Note that particle number conservation is generally not assumed, which 
allows for carrier absorption and emission by the medium (see examples below). 
To derive a quasi-steady state solution for general /, an arbitrary number M + 1 of moments 
are defined: 

^'^^/(x,v) (2) 

d%vfi^,v) (3) 

d'^vvkVif{^,\) (4) 

etc., where h is Planck's constant (spin degeneracy can be included later). Equations for M 
moments follow from multiplication of Eq. ([1]) by 1, f^, f^f/, and velocity integration: 

d^n + V = P(") (5) 
dti + V -n-na = P(^) (6) 

etc. {k,l = l,...,d). The highest order moment and the terms on the right hand side, 
pin) ^ {jn/hY j d'^vCifo - /), P(^') = {m/hYJ d'^vvCifo - /), etc., are functionals of the 
unknown distribution /. In order to close the M moment equations, we will calculate / from 
entropy production rate minimization, subject to constraints of M fixed moments, given by 
Eqs. ([2]), ([3]), etc.. The distribution and derived quantities will thus depend on them. 

Determination of f - It is important to consider the total entropy production rate of the 
whole isolated system, which consists of two contributions associated with the independent 
electron gas and with the medium acting as an equilibrium bath. The first part, sp, can 
be derived from the entropy density ^^^(x) = ~k {^Y J d^v{f In f + (! — /) ln(l — /)} by 
differentiation of sp with respect to time, replacement of dtf with the help of Eq. ([1]), partial 
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f -integration, and finally writing the result as an entropy balance equation, dtSF + V ■ q^^ = 
sf, with entropy current density qs^. This gives sp = k J df^v In {jzj^ ^{f ~^ /o)- 
The second part, sb = W/T, associated with the medium, is obtained from the heat power 
density W = {fY J d'^v (H - /x) C{f - /o). The relation H - fi = A;Tln(l//o - 1) implies 
Sb = k {^Y I d'''v\rL '^{f ~ fo)- The total entropy production rate, sb + sp, is 

The distribution function / is then determined by the variational principle Ss/6f = subject 
to the constraints Eqs. ([2]), Q, ... of fixed moments. This is the main result, and provides 
a closure for the generalized hydrodynamics of the independent electron gas far from local 
equilibrium and for an arbitrary number of moments. 

Two-moment closure - For illustration, consider M = 2 moments and assume C{fo — f) = 
'"(^)(/o ~ /) with f-dependent relaxation rate r. An example beyond this relaxation-time 
approximation will an elastic scatterer in the last example. Optimization of Eq. ([7j) with 
constraints ([2]) and ([3]) leads to 

with Lagrange parameters A*^"'-' and A*^-'^ It is readily checked that the optimization problem 
is convex, which ensures that the solution is a constraint minimum. Solving Eq. ([8]) for / 
and elimination of the Lagrange parameters provides /(v; n,j) and then P^''^\ P^^\ and 11^;. 



Weak Nonequilihrium - Expansion with respect to 5/ = / — /o gives 5f = /o(l — /o)('^^"'' + 
XU) . v)/2r. Elimination of A^") and A^^) with the help of Eqs. ([2]) and ([3]) leads to P^'^) = 
r^"'\nQ — n) and P*^-'^ = — r'^-'^j with 

jd'^vr-\v)U{l-h) ^ ' 

jd'^vv^r-\v)Ul-h) ' ^ ' 

For (i = 3, Eq. (fTOl) gives the usual linear response relaxation-time mobility, e/mr^^'> 
The stress tensor, 11^/, stays diagonal and isotropic, and deviates from equilibrium Llo^fcfc by 

_ {n-n,)jd'^vv\-^h{l-h) 
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Note that /o(l — /o) = —^dyfo, and Eqs. fl9|)- f|TT|) simplify in the low temperature limit 
for /i = mvp/2 > 0, because d^fo is a Dirac 5-function at the equilibrium Fermi velocity 
VF = h{no/ndf''^/m (with fii = 2, ^2 = tt, fig = 47r/3). For T 0, Hq^h = SkmovH {d + 2). 
Finite temperature corrections may be considered in terms of 7 = 2kT/mvp, in the same 
way as done by other standard methods for weak nonequilibrium. 



Fermi Sphere, Zero Temperature - For /i > and T ^ 0, Eq. ([8]) can be solved analyt- 
ically for strong deviation from equilibrium. After multiplying Eq. ([8]) by /(I — /), using 
(1 — /o)//o = exp(m(f^ — Vp)/2kT) and the fact that / G [0, 1], one can solve the resulting 
equation with an expansion / = /(o) +0(72). One concludes for 7 ^ that /(°) = 1 

and Z*^"-* = if /^(v) := m(f^ — Vp)/2kT — X^^^r — A*^-'^ ■ v/r is negative and positive, re- 
spectively (note A*-"'-'-' oc I/7). Hence, i]('v) = defines the boundary of the nonequilibrium 
Fermi surface for T — 0. For f -independent r, the Fermi sphere is shifted according to j/n 
and resized according to n/no as one expects, and the stress tensor becomes 

Hui = ^ <5h + ^ , (12) 

in accordance with the usual expression of the hydrodynamic momentum flux tensor. A 
factor 2 in Qd must be included for spin degeneracy; for instance, Eq. ( fT2l) leads then to the 



d = 3 electron gas pressure introduced by Bloch |10l ]. 

For f-dependent r these results change in general far from equilibrium. As an example, con- 
sider r{v) = tqVp/v with constant tq. For Id, one finds P^"-* = noroln ^/(n/no)^ — {j/nvpY 
and r*^-'-' = r^n^jn. The divergence of P*-"-* at j = vpv? jn^ is an artifact due to the specific 
r{y) and occurs when the shifted Fermi 'surface' approaches f = 0. From P^"\n,j) = and 
Eq. ([6]) one can derive the relation between j and the electric field. A simple experimental 
test of the theory would be based on a measurement of the nonlinear current-voltage char- 
acteristics of a nano-wire on an appropriate support material, if the associated wire-support 
interaction in terms of r{v) is known. 

A non-spherical deformation of the Fermi sphere for varying current can occur in > 1. For 
illustration, we have calculated the deformed Fermi circle in c? = 2 for r{v) = tqVf/v. By 
comparing the result (solid thick curve in Fig. [1]) with a Fermi circle for constant r (dashed 
curve), one observes that the former is closer to the equilibrium distribution for low f- values 
and further away for larger values. This reflects the stronger equilibration, induced by the 



entropy principle [a], at low v for this r{v). It is clear from Eq. (jl]), that not only the shift 
but also the deformation of the Fermi sphere gives a contribution to an anisotropy 11^/. 




FIG. 1: Far from equilibrium Fermi circle in 2d for n = uq and j = nQVF/2 in x-direction for 
constant r (dashed) and deformed circle for r oc l/v (solid). The thin line refers to equilibrium. 

Energy Gap, 3d, Finite Temperature - Consider now fi < 0, such that /o ~ 
N exp{—mv'^ /2kT) <^ 1, with = exp{fi/kT), describes a dilute, non-degenerate, equi- 
librium electron gas in an insulator material. The terms —rf and r/o describe trapping 
and emission of thermally equilibrated electrons from electron traps, respectively. The 
non-equilibrium distribution is calculated for a step function-like rate r{v) = ri for v 
<^ vt '■= \/kT/m and r{v) = r2 = 0.1 ri for v ^ vt, which models a mobility edge 
(inset in Fig. [2]). The resulting distribution function / is shown in Fig. [3] for three different 
cases. Again, in f-regions of larger scattering rates r(t>), / is closer to /o because equili- 
bration is stronger. The effective rate r^^\n,j), far from equilibrium, as a function of the 
current density j (with j = je^ in z-direction) is shown for different values of n in Fig. [2l 
The near equilibrium result, r^^\nQ,j = 0) = 0.306 ri can be obtained directly from Eq. 
( ITOl) . For large currents, j ^ uoVt, as one expects r^^^ r^- It is clear that an additional 
consideration of higher order moments will allow to describe hot electrons. 

Lovj frequency admittance of a Id wire - The generalized hydrodynamics framework is 
very useful to describe low-frequency transport in mesoscopic conductors. To illustrate 
this, the admittance of the one-dimensional symmetric wire shown in FigH] is calculated. 
The contacted long end pieces are described by a finite r = Tc, which models in a natural 
way contacts to macroscopic electron reservoirs, where carriers from the wire are absorbed 
{~'^cf) and equilibrated carriers are emitted (fc/o) into the wire. For —L/2 < x < L/2 the 
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FIG. 2: Mean rate r^^V?"! 

as a function of j for r{v) shown in the inset; dotted: n/riQ = 0.6; sohd: 
n/uQ = 1; dashed: n/riQ = 1.4; dashed-dotted: n/riQ = 1.8. Near equihbrium (upper thin hne): 
r^^^ = 0.306ri; lower thin hne: r2/ri =0.1. 




FIG. 3: Distribution as a function of v in current direction (arrow in inset). Solid: equilibrium 
/ = /o; dashed: j = 0, n = 2no; dotted n = no, j = uqVtI'^- Inset: /-contour for n = no, 
j = UoVt/^. Tail regions with lower r{v) (weaker equilibration) show stronger nonequilibrium. 

electrons are supposed to move ballistic, except for a localized elastic scatterer at x = 
with transmission probability T = 1 — TZ. Hence, C{fo — /) = i'6{x)vF{f{—v) — f{v)) in 
-L/2 <x < L/2 with u = TZ/T. 

Since the purpose is to illustrate the principle, the discussion is restricted to linear re- 
sponse (cf. [ll|). Coulomb interaction can be included self-consistently [l2|. We assume 
negligible capacitive coupling between the contacts, small contact impedance between the 
substrate and the contacted wire pieces, and obviously disregard electron-electron interac- 
tion effects like Luttinger liquid formation [13] or thermal conductance decrease jl^ . 

The low-frequency admittance, G = Gq — iuE, for an applied voltage AV exp{—iujt) , 



with AV = Vi — V2, includes the DC conductance, Go, and the emittance, E [l2|, [iSl]. The 
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FIG. 4: Quantum wire connected to two electrodes at temperature T and potentials Vi,2; with 
contact impedance per length, Zc, and scattering rate r = re- The ballistic region (r 0) of length 
L, contains a localized elastic scatterer with transmission probability T. 

linearized hydrodynamic equations (l5])-(j6]) can be written as 

vpnAn + dxi = (13) 
2e 

K] + vpd^An = —d.^U - 2u6{x)j , (14) 

where dx denotes the spatial derivative. An = n — no, and k = (r — iuj)/vF with r = r{vF)- 
By virtue of = /mo /2m, uq was eliminated. The applied voltage is equal to the total 
difference in the electro-chemical potential, V = U — hvFn/2e. The electrical behavior 
of the wire-support contact is modeled by a phenomenological contact impedance Zc with 
—ezcdxj = Vi,2 — U{x) on the two sides. 

In order to derive G, we integrate Eq. ([Hj) from x = —00 to x = 00 and obtain with 
AV = Vi — V2, I{x) = — ej, and Kc = (rc — iuj)/vF- 

J^V = J Kj{x)dx + {v-^)m, (15) 

where the symmetry of the wire was used. The current in the wire in the left contact 
region x < —L/2 must decay, hence j oc exp(/3a;) there, with [3 = + 2e'^Zc/hKc)~^^'^ 
that follows from Eqs. f|T3l) and f|T4|) . After evaluation of the integral in Eq. f|T5|) one 
obtains e^AV/hI{0) = 1 + z/ + e^ZcuJh - iujL/2vF- With l + u = 1/T and 0, the low 
frequency admittance G = /(0)/A\^ becomes G = ^T+icu^T"^ with D = Ae^L/hvF, where 



factors of 2 for two spin states are included. This result is in accordance with Ref. 121] in the 
considered limit case of vanishing capacitance between the contacts. V{x) in the ballistic wire 
regions has the meaning of a quasi Fermi-level, and the local entropy production rate s{x) = 
mr{P + VF^An^) /Tn^e^ is localized in the contact regions. A generalization to ballistic 
electrons in arbitrary geometries and arbitrarily far from equilibrium is straight-forward; it 
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requires higher order (rnultipole) moments in d > 1, similar to a P — N approximation in 
(non- diffusive) radiation 16|. The generahzed hydrodynamic equations with the discussed 
closure may serve as a general footing for simulations of nano-electronics devices in the full 
range between diffusive and ballistic transport. 



Conclusion - It is also straight-forward to extend the method to other than parabolic 
energy-momentum relations and to generalized moments [l^]. For instance, it should be 
possible to treat in a similar way massless Fermions like neutrinos in stars or electric con- 
duction in graphene, if these particles are independent and the particle- medium interaction 
can be modeled by £(/o — /)• 
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